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Abstract 

T-H ■ 

K^ I We consider the potential benefits of searching for super symmetric dark- 

^^ ' matter through its inelastic excitation, via the "scalar current", of low-lying 

~l I collective nuclear states in a detector. If such states live long enough so that 

f— ^ ' the 7 radiation from their decay can be separated from the signal due to 

^— I I nuclear recoil, then background can be dramatically reduced. We show how 

On ' the kinematics of neutralino-nucleus scattering is modified when the nucleus 

^1 I is excited and derive expressions for the form factors associated with exciting 

^ ' collective states. We apply these results to two specific cases: 1) the I'^ = 5/2"'" 

Qh! state at 13 keV in Ge, and 2) the rotational and hence very collective state 

Q^' I'^ = 3/2"*" at 8 keV in ^^^Tm (even though observing the transition down 

D I from that state will be difficult). In both cases we compare the form factors 

rS^ • for inelastic scattering with those for elastic scattering. The inelastic cross 

^ I section is considerably smaller than its elastic counterpart, though perhaps 

not always prohibitively so. 
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A number of groups are trying to detect weakly interacting dark matter, one of the most 
promising candidates for which is the supersymmetric "lightest neutralino". A popular 
approach is to try to observe the scattering of these particles on nuclear targets in low- 
background laboratory experiments. The signature of neutralino-nucleus scattering is the 
low-energy recoil of the nucleus in a detector. Since the scattering rate is expected to be 
tiny, background is the main factor limiting sensitivity, even when low itself. 

Supersymmetric dark matter is reviewed in ref. |Q] . Here we are interested only in the 
nuclear physics aspects of this problem, and in particular in the possibility of detecting 
inelastic scattering, thereby dramatically reducing background. (The nuclear physics of 
dark matter detection is reviewed in ref. [^.) The work was inspired by questions from 
researchers in the field |^, Q. 

Though inelastic scattering of neutralinos has been considered before, notably in ref. 
[P, the focus was on spin-dependent scattering. The authors discussed low- lying excited 
states in stable nuclei with large measured Ml matrix elements; later, ref. [0 reported 



an upper limit of 9.8 x 10^^ counts/kg/day (at 90% CL) for the inelastic excitation of 
the 7/2+ state at 57.6 keV in ^^^I. It has since become clear, however [0], that spin- 
independent scattering will almost always occur with greater probability than its spin- 
dependent counterpart. We therefore focus here on the possibility of excitation by the 
scalar current, where the relevant multipole is E2 instead of Ml. Collective E2 transitions, 
of which there are many, may allow the scalar current to be even more effective. 

Of course there is a price to pay for the extra 7-ray in the signal from inelastic scat- 
tering: the cross section is noticeably smaller than the elastic one. As we explain below, 
this is caused here not so much by the kinematics discussed in ref. |5[ — E2 excitations 
can often be found lower in the spectrum than Ml excitations — or by the factor qR 
that enters higher multipoles, but rather by a considerable reduction in "coherence" from 
elastic scattering, even when collective nuclear states are excited. Collective excitations 
of the nucleus generally involve valence nucleons, of which there are more than the (effec- 
tively) one that participates in spin-dependent scattering, but still far fewer than the A 
that are involved in elastic scattering. Thus, though we gain in some ways by considering 
the scalar current, we will still not obtain cross sections that approach those from elastic 
scattering. We quantify this remark below. 

Let us consider kinematics first. A particle of mass Mx moves with velocity v and 
scatters on a stationary target of mass Ma- After the scattering the target has -Ecxc of 
excitation energy, i.e. its mass is Mj = Ma + -E-exc- The momentum transfer is 

q^ = Ml\v - v\'^ = Mllv"^ + v'^ - 2vv' cos{e)] , (1) 

where is the scattering angle and v' the final velocity of the scattered particle. The 
energy transfer is 

-» 2 

UJ = Mx {V^ - v'^)/2 = Brecon + ^^exc = 7^ + £^cxc • (2) 

The minimum and maximum momentum transfer, and thus also the minimum and maxi- 
mum recoil energy -Brecoii = 9^/2Mj, correspond to cos(^) = ±1. Eliminating v' we obtain 
a quadratic equation for q'^ which gives 



im = /^^ (^1 ± y 1 - -^ j , (3) 

where // = MxMj /{Mx + Mj) (we can neglect the small difference between Ma and 
Mf here) is the reduced mass. Thus, for the inelastic process to occur at all, we must 
have -Ecxc < [i-v^ l2. (Note that /if^/2 is less than the neutralino kinetic energy, since 
[I < Mx-) To obtain the scattering rate of neutralinos with some velocity distribution 
at a fixed momentum transfer q (or recoil energy -Erecoii)) we have to integrate over the 
velocity distribution from minimum velocity 

^^min = t;- H • (4) 

2m q 



At the same time, for inelastic scattering there is an absolute minimum of momentum 
transfer, q = ^/2JIE^. 

Turning to the nuclear matrix elements that govern the cross section, we have from 
eqs. (4.24), (4.25) of ref. (generahzed to transitions from J ^ J' j^ J) 

^ = (27ffj^^^(^) ' (^) 

where the form factor for initial and final states of the same parity^ is 

Ssiq)= E \{J'\\CLiq)\\J)W (6) 

L even 

and 

ChMiq) = J2'^oJLiqri)YL,M{ri) . (7) 

i 

The summation over L is restricted by \J — J'\ < L < J + J' and the lowest allowed L 
generally contributes most. For appropriate values of J and J' this value will correspond 
to the L = 2 quadrupole mode, which also has the advantage of producing collective 
excitations of the nuclear surface; we denote the associated form factor by 5*2 (g). We have 
lumped into the constant cq all the particle physics aspects of the problem except the 
overall scaling G\. In the ratio of inelastic to elastic form factors the constant cq drops 
out. 

To calculate the matrix elements in eq. (0) we have to know something about the 
structure of the initial and final states. The g — > limit of the matrix element in eq. (1^ 
for L = 2 can be measured in the Coulomb excitation or electromagnetic decay of the 
excited state. The rates of these processes are usually expressed in terms of the quantity 

B{E2,J ^ J') = \{J'\\er^Y2\\J)\'' /{2J + 1) . (8) 

Let us first consider the attractive 9/2"*" — > 5/2"*" excitation in Ge. That isomeric excited 
state at 13 keV has a long half-life (2.95//s) and a rather large B{E2) (23 Weisskopf units 
for the 7-decay transition 5/2+ -^ 9/2+). We make one crude but reasonable assumption 
here: that the transition density for the excitation is concentrated at the nuclear surface, 
as if the excited state were a vibration. Then we have 

B{E2, J^J')c^ eVo^^^ang)V(2J + 1) • (9) 

where R is the nuclear radius, po the proton density, and Aang the matrix element of the 
angular factors. 

With the same assumptions we can write the form factor S2{q) for the inelastic neu- 
tralino J ^ J' scattering as 

S2iq) « cl\{f\\h{qr)Y,\\J)\' = eg (^^p,'^' j,{qRf{A,^,f , (10) 



^One might imagine i51-like transitions between low- lying states of opposite parity, but for nuclear- 
structure reasons their strengths are notoriously small. 



where the factor AjZ conies from the additional assumption that the neutron and proton 
densities are proportional. Using the known B{E2) we can rewrite the above as 

S2{q) = cl^i2J + l)h{qR?^^ ■ (11) 

The 5*2 form factor can then be compared to the form factor for elastic scattering, which is 
governed by the operator Cqo = coX^j Jo('7'"i)^oo(^j)- A constant density inside the nuclear 
radius and the relation 

/ joiqrydr = —jiiqR) (12) 

Jo q 

give 

S.iq)=cl{2J+lfA^'jl^. (13) 



The ratio of inelastic to elastic cross sections, S2{q)/Sei{q), from eqs. (|ll| ) and (|l3), is 
independent of the constant cq . 

Figure 1 (the upper panel) shows the elastic and inelastic form factors as a function of 
the recoil energy -Erecoih normalized to the elastic form factor at q = (i.e. -Erecoii = 0). The 
inelastic form factor in fact begins at a finite -Brecoii related to the minimum momentum 
transfer in eq. (^). The largest Ey-ecoH we consider, 140 keV, corresponds to neutralinos of 
mass ~ 60 GeV, (the mass indicated by a recent experiment [^) moving with the galactic 
escape velocity, 650 km/s. For inelastic scattering, eq. (y) restricts the i^recoii to less than 
about 127 keV. At low recoil energies, -Brccoii ^ 30 keV, the inelastic form factor is small 
because the spherical Bessel function J2{x) is proportional to a;^/15 for small x. Even at 
larger recoil energies, however, the inelastic form factor is down from the elastic one by 
a factor of 100 - 1000. Only near the zero of the function ji{x), which corresponds to 
-E'recoii ~ 220 keV in Ge, is the inelastic cross section larger than the elastic one. The 
small inelastic cross section is caused by the absence of the coherence factor A^ (which 
appears divided by Z in eq. ( |lll) only to renormalize the density). The collectivity of 
the E2 transition, which as noted above is restricted to the nuclear surface, cannot fully 
compensate this loss. Thus, while the sharp 7 ray in the signal is undeniably beneficial, the 
expected count rate is substantially smaller than in elastic scattering. To further quantify 
this statement we evaluate the total elastic and inelastic cross sections for neutralinos with 
Mx = 60 GeV and a Maxwellian velocity distribution (v = 220 km/s) terminated at the 
galactic escape velocity (650 km/s). The result for an ideal detector is 

/ —inelastic \ 

T^S^ = 2-8 X 10-^ . (14) 

A real detector will have some threshold in recoil energy below which it is not sensitive. 
The elastic form factor is largest at low recoil while the inelastic form factor is completely 
negligible there; excluding events with energies below the lower limit will therefore increase 
the ratio above. In a detector with a 10 keV threshold, the ratio is 
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= 5.7 X 10"^ , (15) 
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Figure 1: The quantities S = cr((7)/(7eiastic(9 = 0) for elastic (dotted lines) and inelastic 
(full lines) neutralino scattering. The upper panel (a) is for '^^Ge and the lower panel (b) 
for ^^^Tm. The dashed line in (b) is the inelastic S evaluated with eq. ([ll]), which is less 
accurate than eq. ([l9|). 



still a rather small number. 

Are there circumstances in which the reduction is not so dramatic and an experiment 
more desirable? For this to be the case, there must exist a low-lying (not much more than 
20 keV) excited state with a very collective E2 transition. This state must live sufficiently 
long so that its deexcitation can be separated in time from the signal caused by the recoil 
kinetic energy. Finally, to eliminate the need for isotope enrichment, the target nucleus 
should be the only stable isotope of the element it represents. 

A quick search of the Table of Isotopes]^] reveal that these conditions are not so easy 
to fulfill. In fact, we found only one nucleus, ^^^Tm, that comes close. Its rotational 3/2"'" 
state at 8.4 keV has a half-hfe of 4.1 ns and very collective B{E2; 3/2"'' -^ 1/2+^.) of 226 
Weisskopf units. Detecting inelastic scattering to this state will be difficult; its excitation 
energy is too low and its half-life too short. Nevertheless, we evaluated the corresponding 
form factor to see what kind of count rates we could expect. 

In nuclei with permanent deformation the B{E2) values are related to the expecta- 
tion value of r^l20 in the intrinsic frame of the nucleus, which in turn follows from the 



deformation parameter /3: 

(r^l1o^n = ^/3(l + ^/J/3 + ...) . (16) 

We can also write the intrinsic-frame expectation value of the operator C20 in eq. (|^) in 
terms of (3: 

(4^(9)) = ^/? (j2{qRo) + ^^[qRoJiiqRo) " j2{qRo)] /? + • • •) • (17) 

To relate the inelastic form factor to the B{E2) value, we use the expressions for rotational 
states: 

B{E2;J,K -^ J',K) = {r'^Y^f')'^ {J K2Q\J' K)"^ , (18) 

and 

S2{q;J,K^ J',K) = {CiT{q)f{2J+l){JK20\J'Kf . (19) 

(The quantum number K, the angular-momentum projection on the nuclear symmetry 
axis, is 1/2 for ^^^Tm.) To leading order in /3, these relations give the same result as eq. 
(PT|). The terms of order /3^ supply about a 10% correction in ^^^Tm, which has P ~ 0.3. 
The form factors for ^^^Tm appear in the lower panel of Fig. 1. The maximum recoil 
energy for a 60 GeV neutralino with the galactic-escape velocity is now 112 keV. The 
inelastic form factor, as expected, is not as suppressed compared to its elastic counterpart 
as in Ge; the factor is less than 100 in the broad maximum of the inelastic form factor 
at ~ 30 keV recoil energy. The ratio of the total cross sections integrated from the lowest 
possible momentum transfer is now 

/-j.inelastic\ 

\^elastic) = 1-5 X 10-' , (20) 

and increases to 

/^inclasticv 

/^elastic/ = 5-9 X 10-^ (21) 

\" / from 10 keV 

when integrated from a 10-keV threshold. To relate these results to those in Ge, one must 
recall that the normalizing factor, Sei{q = 0), scales like A^, i.e. it is larger for ^^^Tm than 
for ^^Ge by (169/73)^. With a 10 keV threshold, the integrated inelastic cross section per 
kg of material in ^^^Tm is therefore suppressed with respect to the elastic cross section in 
^■^Ge by less than 100, a number that may not be so intimidating. 

In conclusion, we have examined neutralino inelastic scattering to collective states 
with large B{E2) values. We have shown how to evaluate the form factors and presented 
examples. While the search for inelastic neutralino scattering offers an opportunity to 
suppress most background, it also leads to a considerable reduction of the expected signal. 
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